Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN WF$t style file v2.2) 



On the relation between electron temperatures in the + 
and ++ zones in high-metallicity H II regions 



o ; 
o ■ 

(N ■ 

a ■ 



Leonid S. Pilyugin 1 , Jose M. Vflchez 2 and Trinh X.Thuan 3 

1 Main Astronomical Observatory of National Academy of Sciences of Ukraine, 27 Zabolotnogo str., 03680 Kiev, Ukraine (pilyugin@mao.kiev.ua) 
2 Instituto de Astrofisica de Andalucia, CSIC, Apdo, 3004, 18080 Granada, Spain (jvm@iaa.es) 

3 Astronomy Department, University of Virginia, P.O.Box 3818, University Station, Charlottesville, VA 22903 (txt@virginia.edu) 



Accepted 2006 May 25. Received 2006 May 25; in original form 2006 April 13 



OS 



> 

in 
o> 
\o 
in 
o 
\o 
o 

6 

u 



ABSTRACT 

We suggest a new way to establish the relation between the electron temperature t% 
within the [O in] zone and the electron temperature ti within the [O Ii] zone in high- 
metallicity (12+log(0/H) > 8.25) Hn regions. The t 2 - t3 diagram is constructed by 
applying our method to a sample of 372 Hn regions. We find that the correlation 
between t 2 and t3 is tight and can be approximated by a linear expression. The new t% 
- t% relation can be used to determine t% and accurate abundances in high-metallicity 
Hn regions with a measured t3. It can also be used in conjunction with the ff relation 
for the determination of t3 and t 2 and oxygen abundances in high-metallicity Hn 
regions where the [OIII]A4363 auroral line is not detected. The derived t 2 - £3 relation 
is independent of photoionization models of Hn regions. 
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1 INTRODUCTION 

Accurate abundances in H 11 regions can be derived via the 
classic T c method, T e being the electron temperature of 
the Hn region. In this method, the reliability of the equa- 
tions for determining the ++ /H + and + /H + ionic oxy- 
gen abundances dep ends to a great par t on th e reliability 
of the atomic data. iPilvuein fc Thuan I i2005l) have com- 
pared the ionic oxygen abundances derived with a recent 
set of equations using the latest atomic data (Ilzotov et alJ 
( 2006) ), with thos e dete rmined from an earlier set of equa- 
tions iPaeel et alJ 1)19921) ) and have found close agreement. 
Thus, the equations for the determination of the ionic oxy- 
gen abundances in Hn regions appear to be robust, in the 
sense that, for a given set of measured electron temperatures 
t3 and t2, relations determined by different authors result in 
close oxygen abundances in H 11 regions. 

The electron temperatures £3 within the [O ill] zone and 
ti within the [0 11] zone are determined from diagnostic line 
ratios when they are available. The intensity ratio of the 
[O IIl]A4959 + A5007 nebular line to the [O IIl]A4363 auroral 
line is used to determine t3. As for ta, the intensity ratio of 
the [O ll]A3727 nebular line to the [0 11] A7320+A7330 auroral 
line is used. Based on a sample of H 11 regions where all the 
necessary oxygen lines are detected. iKennicutt et all J2003T) 
found that the ta temperatures show a large scatter and that 
they are nearly uncorrelated with t3 temperatures. They 
noted that the measurements of the faint [0 11] A7320+A7330 



auroral l ines may contain large rand om errors. Furthermore, 
iRubin I Jl986f ): iTsamis et all (120031) have pointed out that 
the [On1A3727/[O ii1A7320 + A7330 r atio may be affected 
by recombination. Ilzotov et alJ ll2006h have derived ta and 
t3 from the above ratios for a sample of Hn regions. They 
concluded that the two temperatures generally follow the 
t2 - t3 relation obtained from photoionization models (e.g. 
Stasinska 1982, 1990), but that the scatter of the data points 
is very large. The large scatter was attributed to large flux 
errors of the weak [O n]A7320 + A7330 emission lines. 

The intensity ratio of the [Nll]A6548 + A6584 nebu- 
lar line to the [Nll]A5755 auroral line is also used to de- 
termine t2. It seems to give more reliable values. How- 
ever, the contribution of recombination to the excitation 
of the [Nll]A5755 line may affect tempera tures derived 
from the [Nil] neb ular to auroral line ratios iRubin Ill986l : 
ITsamis et all 1200 j) . The ([N n]A6548 + A6584)/[N n]A5755 
ratio is measured in only a few H 11 regions with a detected 
([Om]A4959 + A5007)/[O m]A4363 line ratio. 

When only a single temperature measurement is avail- 
able, a ta - t3 relation based on grids of H n region mod- 
els is usually used. Several versio ns of such a £2 - £3 
relation have been proposed (e.g. JCamrjbell^^l 

iDa 



Paeel et al]ll992t Ilzotov et al.lll997t iDeharveng et alJl200C 
Oev fc Shields! l2000ft . but the agreement between them is 



not very good. The available measurements do not provide 
an undisputable evidences in favour of any out of suggested 
relations, i.e. the choice of the relation is in fact arbitrary. 
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Thus, the £2 - £3 relation seems to be the weakest link of 
the classic T c method. We will address this problem here. 

The main goal of this paper is to derive a t2 - t3 relation 
which is model-independent. This allows one to relax the ar- 
bitrarity in the choice of the £2 - £3 relation. The basic idea 
is the following. Usually, the set of equations used for the de- 
termination of the ionic ++ /H + and + /H + abundances 
in H 11 regions is applied to the whole nebula. If we apply 
the equation for ++ /H + only to the ++ zone, then this 
would yield, not the ionic ++ /H + abundance, but the total 
O/H oxygen abundance instead. There is no reason to sus- 
pect that the oxygen abundances in the ++ and + zones 
differ. We therefore require that the equation for ++ /H + 
applied to the ++ zone and the one for + /H + applied to 
the + zone result in exactly the same value of the oxygen 
abundance. This condition allows us to derive a relation be- 
tween t2 and t3. We can apply the equation for ++ /H + to 
the ++ zone alone, only if we know the contribution of the 
++ zone to the measured Hg flux of the Hn region. We 
will b e using the ff relation <PilvugirJl2005l IPilvugin et al.l 
2006) for this purpose. 

The observational sample to be used is described in Sec- 
tion 2. The strategy for determining the t2 - t3 relation is 
discussed in Section 3. A model-independent t2 - t3 relation 
is derived in Section 4. We discuss the results in Section 5, 
and summarize our conclusions in Section 6. 

We will be using the following notations through- 
out the paper: R 2 = I[oii]a3727+a3729/Ih /3 , R3 = 

I[OIII]A4959+A5007/ 1 H /3 , R = I[OIII] A4363/Ih 3 , R23 = R2 + R3- 

With these definitions, the excitation parameter P can be 
expressed as: P = R3/(R2+Rs). 



2 OBSERVATIONAL DATA 

A large sample of high-precision measurements of H II re- 
gions is at the bas e of the present investigati on. We include 
first the sample of lPilvugin fc Thuan I J2005I) who have car- 
ried out an extensive search of the literature to compile a 
list of more than 700 individual spectra of H II regions in ir- 
regular and spiral galaxies, with the requirement that they 
all possess a detected [O in] A4363 emission line. While we 
have tried to include as many sources as possible, we do 
not claim our search to be exhaustive. Since the majority 
of extragalactic Hn r egions are in the low -density regime 
JZaritskv et alJ Il994t iBresolin et "all l2005h . only such Hn 
regions will be considered here. Of the objects with the 
electron density-sensitive ratio r n = [S n] A6716/[Sll] A6731 
available, those with r n < 1.3 were excl uded. 

W e have also included the sample of llzotov et al.l J2004L 
2006) who have extracted from the Data Release 3 of the 
Sloan Digital Sky Survey (SDSS) around 4500 spectra of 
H II regions with an [O ill] A4363 emission line detected at a 
level better than la, and carefully measured the line inten- 
sities in each spectrum. The SDSS Hn regions are also in 
a low-density regime. Yuri Izotov and Natalia Guseva have 
kindly provided us with the total list of t heir measurements, 
as on ly part of these have been published lllzotov et al. 2004, 
2006). In total, our sample consists of around 5200 Hn re- 
gion spectra. 

Followi nglPilvugin fc Thuan I J2005I) , we can use the the 
ff relation JpilvuginetHdT 2006T) to select out the H II re- 



gions with the highest quality measurements. We can define 
a "discrepancy index", equal to the difference between the 
logarithm of the observed flux R obs in the [Om]A4363 line 
and the logarithm of the flux R cal of that line derived from 
the strong [Oll]A3727, [O m]AA4959,5007 lines using the ff 
relation: 



D a =lo gj R° 



log R c 



(1) 



It is well known that the relation between the oxy- 
gen abundance and the strong oxygen line intensities 
is double- valued, with two distincts parts, traditionally 
known as the upper high-metallicity branch and the lower 
low-me tallicity branch of t he R 23 - O/H diagram. Fol- 
lowing IPilvugin fc Thuan I <|2005f) . we adopt the value of 
12+log(0/H) = 8.25 as the boundary between the upper 
branch and the transition zone. The exact boundary is dif- 
ficult to establish, but we have chosen that value because 
we will be using later the ff relation which is applicable at 
metallicities above 12+log(0/H) ~ 8.25. We have extracted 
from our total sample a subsample of high-metallicity H 11 re- 
gions with an oxygen abundance 12+log(0/H) > 8.25. The 
o xygen abundances h ave been calculated with the equations 
of llzotov et alJ i200(Jl for the T e m ethod, and with the t2 - 
fa relation of lCamabelTet alJ lll986F) . Fig.Qshows the cumu- 
lative number of individual Hn region measurements with 
the absolute value of the discrepancy index Dff less than a 
certain value. Solid and open circles show data for Hll re- 
gions with respectively positive and negative values of the 
discrepancy index. It is seen that the agreement between the 
cumulative numbers of measurements with positive and neg- 
ative values of the discrepancy index Dff is good for values 
of Dff less than 0.05, but gets increasingly worse at higher 
values. This is the result of a selection effect. Indeed, a large 
positive value of Dff means that the observed R obs is signif- 
icantly overestimated. This leads to a substantial overesti- 
mate of the electron temperature and an underestimate of 
the oxygen abundance, bringing it lower than our metallicity 
cut-off of 12+log(0/H) = 8.25. This effect thus systemat- 
ically depletes the number of high-metallicity Hll regions 
with a large positive Dff. To free ourselves from this system- 
atic effect, we select only measurements with the absolute 
value of Dff less than 0.05. This sample contains 372 data 
points and will be referred to as the standard sample. The 
H 11 regions in it have line intensity errors that are random. 



3 A STRATEGY FOR THE DETERMINATION 
OF THE t 2 - t 3 RELATION 

3.1 The usual version of the T e method 

For abundance determination in an H 11 region, a two-zone 
model for its tem perature structure is usually adopted. 
Ilzotov et alJ i2006l) have recently published a set of equa- 
tions for the determination of the oxygen abundance in H 11 
regions in the context of such a two-zone model. Accord- 
ing to those authors, the electron temperature £3 within the 
[O ill] zone, in units of 10 4 K, is given by the following equa- 
tion 

(2) 



t 3 



log(R 3 /R)-\ogC T ' 



The quantity Ct is defined by: 
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ta = 0.243 + 1.031 1 3 - 0.184 1| 



(12) 



0.00 0.05 0.10 0.15 0.20 0.25 
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Figure 1. Cumulative number of individual measurements of H II 
regions with the absolute value of the discrepancy index Dff less 
than a given value. The solid circles show data for H II regions with 
positive values of the discrepancy index, and the open circles for 
H II regions with negative values of the discrepancy index. 



Based on H n r egion model calculatio ns by 
IStasiriska fc Schaererl Jl997fl . iDeharveng et alJ (120001 
have suggested the following relation 

t 2 =0.775 t 3 +0.281. (13) 

lOev fc Shields! (12000) have found that the Campbell et al. 
relation is reasonable for *3 > 1.0. However at lower temper- 
atures, the models are more consistent with an isothermal 
nebula. They consequently adopted the formulation, 



0.7* 3 + 0.3, 

*3, 



t 3 > 1.0 
t 3 < 1.0. 



(14) 



Using the system of equations Eqs.Q-i-© and one of 
the t2 - t3 relations chosen among Eas. llUI - M14l '). the oxy- 
gen abundance (0/H)t o in an Hn region can then be deter- 
mined. 



C T = (8.44 - 1.09 t z + 0.5 t% - 0.08 t%) v, 

where 

_ 1 + 0.0004 x z 
" ~~ 1 + 0.044a;3 ' 



and 

X3 = 10 _4 n e ^ 



-1/2 



(3) 
(4) 
(5) 



As for the ionic oxygen abundances, they are given by 
the following equations: 



12 + log(0 ++ /tf+) =log(/ [om] 



A4959+A5007/-Th / j) + 
1 251 

6.200 + - 0.551og* 3 -0.014*3, (6) 

t 3 



and 



5.961 + 



12 + \0g(0 + /H + ) = log(/ [0 //]A3727+A3729/^H f ,) + 

-4-^- -0.401ogt 2 - 0.034 i 2 +log(l + 1.35K2), (7) 

*2 



where 

,„-4 ,-1/2 

X2 = 10 n e t 9 



(8) 



Here n e is the electron density in cm" . 

The total oxygen abundances are then derived from the 
following equation: 

H H+ + H+ ' ( ' 

The electron temperature t 2 of the [O n] zone is usually 
determined from an equation which relates t 2 to * 3 , derived 
by fitting Hn region models. Several versions of this t 2 - 
£3 relation have been proposed. A widely used relation is 
the one bv lCampbell et alJ ljl98rt Csee also lOarnettJ ()l99^ 1 
based on the Hn region models of Stasiriskaj 1119821) : 

t 2 = 0.7* 3 + 0.3. (10) 

Another relation has been propose d by [ Pag el el^jd.1 iflC 
also based on H 11 region models of IStas iriska Jl99fi) 1 

1=0.5(1+0.8). 

t 2 13 



(11) 



Izotov et al 1 lll997l) . fitting also the Hn region models of 
Stasiriskal <199Q) . have proposed the following expression 



3.2 An alternative version of the T c method 

Examination of Eq.© shows that the fluxes used to calcu- 
late R3 and in the H/3 line originate from different volumes 
of the nebula. The H/3 flux comes from the whole Hn region 
while the R3 fluxes come from the volume where the oxygen 
is in the ++ stage. If Eq.JHJ is applied to the ++ zone 
only, then it gives the total Q instead of the ionic oxy- 
gen abundance as in the traditional approach. In that case, 
Eq.@ takes the form 



12 + \og{0/H) 

+6.200 



log(%5//I]A4959+A5007/(™ X Ih^)) 

1.251 



t 3 



0.55 log t 3 -0.014*3, (15) 



where w is the fraction of the Hp flux in the ++ zone. With 
the adopted notations, Eg. (1151 can be rewritten as 



12 + log(0/fT) 
1.251 



ta 



= log(R 3 /w) +6.200 
0.55 log ta -0.014*3. 



(16) 



Similarly, Eq.JTJ can be rewritten as 



12 + \og{0/H) 



I°g(^[OiT] A3727+A3729/ [(1 

1.676 



+5.961 + 



*2 



w) X I H$ ] 

- 0.40 log t 2 - 0.034*2 

+ log(l + 1.35x2). (17) 



12 + \og{0/H) = Iog[Ba/(l - w) 



5.961 



1.676 

ta 



-0.401og* 2 - 0.034*2 + log(l + 1.35k 2 ). 



(18) 



There is no reason to expect the oxygen abundance in the 
++ zone to differ from that in the + zone within the same 
Hll region. In other words, the oxygen abundance derived 
from Ea. l|16|l must be equal to the one derived from Ea. l|18|l . 
Equating the right-hand sides of these two equations result 
in 



log(-R 3 /iu) + 6.200 + 



1.251 

T 3 



0.55 log t 3 -0.014*3 



1 676 

log[R 2 /(l-w)] + 5.961 + 0.401og* 2 

*2 

-0.034*2 + log(l + 1.35x2). (19) 
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Figure 2. The R vs. R3 diagram. The circles show individual 
H II region measurements from the standard sample. The relation 
corresponding to Eq. I2UI for the value of the excitation parameter 
P = 1.0 is shown by the solid line. The dotted lines show from 
right to left the relations corresponding to that equation for P = 

0. 8, 0.6, 0.4 and 0.2. The dashed lines show from top to bottom 
the positions of H II regions with a given value of the electron 
temperature t 3 = 1.2, 1.1, 1.0, 0.9, 0.8 and 0.7, in units of 10 4 K. 

Eg. 1191 can be solved for t 2 if the values of t3, n e and w 
are known. The majorit y of extragalactic H 11 regions are in 
the lo w- density regime jZaritskv et al.lfl994l : iBresolin et alJ 
2005). Therefore we adopt an electron density n e = 100 cm -3 
for all the Hn regions in our sample. Thus, Eo. 1191 allows 
to derive the electron temperature t 2 within the [O n] zone 
in a Hn region if the electron temperature t 3 within the 
[O ill] zone and the fraction w of H,3 flux in the ++ zone 
are known. 

3.3 Estimation of the parameter w 

The fraction w of flux in the ++ zone can be esti- 
mated in the following way. We have investigated recently 
the relationship between the observed auror al and nebular 
oxygen line fluxes in spectra of H 11 regions JPilvugirJl200i 
iPilvuein et aU2 006l . We have found a relation (called the ff 
relation) that is metallicity-dependent at low metallicities, 
but independent of metallicity above 12+logO/H ~ 8.25, 

1. e. there is one-to-one correspondence between the auroral 
and nebular oxygen line fluxes in spectra of high-metallicity 
Hn regions. Using a compilation of recent high-precision 
measurements of oxygen lines fluxes in high-metal licity H 11 
regio ns, the following ff relation was derived iPilvuein et alJ 
2006) 

log 7? = -4.151 -3.118 log P + 2.958 log R 3 

- 0.680 (log P) 2 . (20) 

Fig.|5|shows the relation corresponding to Eg. 1201 in the R3 
- R diagram for various values of the excitation parameter P. 
The solid line, called hereafter as the basic line, corresponds 
to P = 1.0, while the dotted lines correspond, from right to 
left, to P = 0.8, 0.6, 0.4, 0.2. The dashed lines show the loci 
of H n regions for several values of the electron temperature 
t 3 . From top to bottom, t 3 = 1.2, 1.1, 1.0, 0.9, 0.8, 0.7 in 



units of 10 4 K. The circles show individual measurements of 
Hn regions from the present sample. 

Let us consider a sequence of H 11 regions with the same 
value of t 3 . The fraction of radiation of the Hn region in 
the ++ zone increases along this sequence from low to high 
values of P, and reaches its maximum value, w = 1, at P = 
1, i.e. at the intersection point of a dashed line with a given 
electron temperature with the basic solid line . Then, the 
fraction of radiation in the ++ zone can be estimated as 

w = m hs /R% =1 - (21) 

Such an estimate is not airtight since not only w can change 
along a sequence of Hn regions with a fixed value of t 3 , but 
other characteristics of Hn regions, such as, for example, 
their oxygen abundance, can do so as well. Therefore, the 
value of t 2 derived from Eq. 1191 and Eq. 1211 may carry some 
systematic error E t . That error should be small for high- 
excitation Hn regions (they are located close to the basic 
line) but increase with decreasing P. In the following, we 
discuss how to estimate that error and correct t 2 for it. 



4 THE t 2 -t 3 RELATION 

The top panel in Fig.|3]shows the value of t 2 as derived from 
Eg. 1191 and Eg. 1211 for the objects in our standard subsam- 
ple of high-metallicity Hn regions. As noted above, t 2 has 
a systematic error E t that depends on the value of the exci- 
tation parameter P. To illustrate the point, we show in the 
bottom panel of Fig.[3]high quality data (those with |Dff | < 
0.01) for Hn regions with 0.9 > P > 0.8 (filled circles) and 
for those with 0.5 > P > 0.4 (open circles). Inspection of the 
bottom panel in Fig.|S|shows that indeed the t 2 - t 3 relations 
are slightly different for Hn regions with different values of 
P: at a given value of t 3 , t 2 = t 2 + E t is slightly smaller, on 
average, for low than for high excitation Hn regions. The 
true t 2 - t 3 relation can be found by extrapolation of the t 2 
- t 3 relation to P = 1. 

The bottom panel in Fig. [3] shows that the systematic 
error E t is small. Therefore, we will approximate it by a 
linear expression, i.e. the true t 2 and the derived t 2 values 
are related by 

t 2 =t* 2 -E t = t* 2 -c(l-P). (22) 

Examination of Fig. |3] shows also that the i 2 - t 3 relation 
can be parametrized as 

t 2 = at 3 + b. (23) 

The coefficients a and 6 in Bo, (1231 and the coefficient c in 
Eo. 1221 can be found simultaneosly by fitting the data in 
the top panel of Fig. [3J by an expression of the type 

t* 2 = at 3 +b + c (1 - P). (24) 

Carrying out a standard least-squares fit to the data, we 
obtain the following values of the coefficients; a = 0.716 ± 
0.012, b = 0.264 ± 0.013, c = -0.042 ± 0.007. Then, the true 
t 2 - t 3 relation is 

t 2 = 0.716(±0.012) t 3 + 0.264(±0.013). (25) 

The value of w given by Eo. 1211 is the fraction of ra- 
diation in the ++ zone if the condition t 2 =t 3 holds. This 
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Figure 4. The difference AtJ, denned in section 4, as a func- 
tion of excitation parameter (top panel), electron temperature t3 
(middle panel), and discrepancy index Dg (bottom panel) for our 
standard sample of H II regions. 



Figure 3. Top panel. The electron temperatures t| versus t3 for 
our standard sample of Hn regions. Bottom panel. The same as 
in the top panel, except that only high-quality data with |Dfj| < 
0.01 are included. Filled circles show Hn regions with 0.9 > P > 
0.8, and open circles those with 0.5 > P > 0.4. 



is not however the case. Since the Hg emission coefficient, 
E2 2, is dependent on the electron temperature l)Allerlll984h 



E° 2 = 1.387 T 



10 



-0.0424/t 



(26) 



then the contribution w c of the [OIII] zone to the total flux of 
the nebula in the H,3 line is given by the following expression 



#4,2 (*3) 



E° (t 3 ) w + E9 a {ta) (1-w)' 



(27) 



We have recomputed the values of t^j using w c instead of w. 
We found the differences between the t| values derived with 
w and Wc to be less than 0.01. Therefore, we will neglect 
this small correction, and will use the value of w given by 
Ea.pTjl. 

We now check whether the linear forms of the analytical 
expressions adopted for the t2 - t3 and t| - t2 relations are 
justified. The top panel of Fig. 2] shows the difference At 2 
between the value of t| derived from Eg. 1191 and Eg. 1211 
and the one derived from Eg. 11241 as a function of the exci- 
tation parameter P. The points are individual measurements 
of H 11 regions from the standard subsample, and the line is 



the linear best fit to those data obtained through the least- 
squares method. The middle panel in Fig.|l]shows the same 
difference At2 as a function of the electron temperature t3. 
Examination of the top and middle panels of Fig. 0] shows 
that At2 does not coorelate either with P or t3, justifying 
our adoption of linear forms for the t2 - t3 and t| - t2 rela- 
tions. 

The bottom panel in Fig. 2] shows Atjjj as a function of 
the discrepancy index Dg. There is an anticorrelation of the 
two quantities. Since the discrepancy index appears to be an 
indicator of the error in the auroral line R measurements, 
this suggests that the t2 - t3 correlation is rather tight, 
and that any scatter in this correlation is caused mainly by 
uncertainties in the measurements. 

The low uncertainty in the values of the coefficients in 
Eo. 1251 may be a consequence of the large number of points 
(372) used. The following consideration can tell us some- 
thing about the true accuracy of the derived t2 - t3 relation. 
We have extracted from our standard sample a subsample 
of 86 high-quality measurements of Hn regions with |Djt| 
< 0.01. The derived t| - t3 diagram for this subsample is 
shown in the top panel of Fig. |S] A fit to those data gives 

t* 2 = 0.729 (±0.006) i 3 + 0.257 (±0.007) 

-0.058 (±0.003) (1 - P) (28) 

The corresponding t2 - t3 relation is then 

t 2 = 0.729 (±0.006) i 3 ± 0.257 (±0.007). (29) 
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Figure 5. Open circles show the tj - t3 diagram for the sub- 
sample of Hn regions with Dff < 0.01 (top panel), and for the 
standard sample (bottom panel). 



This relation is shown by the solid line in the top panel of 
Fig. For comparison, the bottom panel of Fig. |K| shows 
the t^ - t3 diagram and the t2 - t3 relation for the whole 
standard sample. 

Fig.^shows the differences AtJ between the values of t2 
derived from Eg. 1191 and Eo. l2H and the one derived from 
Eg. 12811 as a function of P (top panel), t3 (middle panel) and 
Dff (bottom panel). The points are individual Hn regions 
with 0.01 > Dff > -0.01, the lines are linear least-square fits 
to those data. Comparison of Fig.^Jand Fig.HJshows clearly 
that the scatter in t2 at a fixed value of t3 decreases in the 
subsample of H n regions with more precise measurements. 
This confirms our conclusion that the t2 - t3 correlation is 
rather tight, and that the scatter in it is caused mainly by 
measurement uncertainties. 

Examination of Eo. 12511 and Eg. 1291 shows that the t2 
- t3 relations derived from the two different samples of H n 
regions are very similar and agree within the formal uncer- 
tainties. Thus the derived t2 - t3 relation is rather robust. 
In the following, we will adopt as the t2 - t3 relation 



ti = 0.72 1 3 +0.26 



(30) 



We note that it has been generally accepted that there 
is a one-to-one correspondence between t2 and t3, i.e. that 
the t2 - t3 relation does not depend on an additional pa- 
rameter. If there is a dependence of the t2 - t3 relation on 
the excitation parameter (from general considerations, this 
cannot be excluded), that dependence cannot be revealed 
by our approach and can influence the derived relation. 




Figure 6. The difference Atj, defined in section 4, as a function 
of excitation parameter P (top panel), electron temperature t3 
(middle panel), and discrepancy index Dff (bottom panel) for the 
subsample of Hn regions with |Dg | < 0.01. 



In summary, there is a tight correlation between the 
electron temperature tz within the [O ill] zone and the elec- 
tron temperature ti within the [O n] zone in high-metallicity 
Hn regions. This correlation can be well approximated by a 
linear expression and its form is rather robust. 



5 DISCUSSION 

We now compare the t2 - t3 relation obtained here with 
those obtained by other authors. Fig. |7] s hows our t2 - t 3 
relation tog e ther w i th those from |Cajm3 belLet_alJ lj986j); 
iPaeel et alJ (I1992T): lizqtov et alJ (Il997ft : iDeharveng et all 
fcOOdklOev fc Shields! (l200(tl . Since our t 2 - t 3 relation is 
derived for cool high-metallicity Hn regions then the high- 
temperature low-metallicity part of relation is not consid- 
ered here. Fig. |7| sh ows that our t2 - t3 r elation is most 
similar to the one bv lCampbell et alJ ([,1986). 

We consider next how the obtained t2 - t3 relation may 
affect the derived oxyge n abundances in Hn re gions. Since 
the t2 - t3 relation of ICampbell et al.l lll986l) has found 
wide acceptance and use in abundance determinations in 
Hn regi ons, we will compare oxygen abundances derived 
with the lCampbell et al.l dl986fl relation and ours. For this 
purpose, we use the dat a for H u regions in th e spiral galaxy 
NGC 4 254 obtained bvlMcCall et all Jl985tl : IShields et all 
Jl99ll) : iHenrv et alJ lll994h . Since the auroral [Om]A4363 
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Figure 7. Comparison of the t2 — t3 relation derived here with 
those derived by other investigators. 
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Figure 8. The radial distribution of oxygen abundances across 
the disk of the spiral galaxy NGC 4254. The filled circles show 
(0/H)ff abundances based on the — t3 relation derived in the 
present study, and the solid line is the linear least-squares best fit 
to those data. The open circles show (Q/H)ff a bundances based 
on the t2 - t3 relation of Ca mpbell et all Jl98fit) . and the dashed 
line is the linear least-squares best fit to those data. 

line was not detected in those Hn regions, 0/H)g abun- 
dances based on the line fluxes R cal have been determined 
from the ff relation, Eq. (l20H . for every Hn region, following 
IPilvugin et all (l2006l) . The (0/H) ff abundances in the H II 
re gions of NGC 4254 det ermined using the %2 - t3 relation 
of ICampbell et ail |l986) are shown as a function of galac- 
tocentric distance in Fig.[S]by open circles. The dashed line 
is the linear least-squares best fit to those data: 

12 + \og(0/H) = 8.74 (±0.03) - 0.036 (±0.004) x R g . (31) 

We have adopted a distance of 16.14 Mpc for NGC 4254 
flPilvugin et al.ll2004h . The filled circles are (0/H) ff abun- 
dances determined with our t2 - t3 relation. The solid line 
is the linear least-squares best fit to those data: 



12 ± \og(0/H) = 8.86 (±0.03) - 0.045 (±0.004) x R g . (32) 

Fig. |H] shows that the abundances derived with our t2 - t3 
relation are slightly higher (up to ~ 0.1 de x ) tha n those 
derived with the relation by ICampbell et all jl986T) . Com- 
parison of Eg. 1311 and Eg. H32H leads to the same conclusion. 



6 CONCLUSIONS 

We suggest a new way to establish the relation between the 
electron temperature £3 within the [O ill] zone and the elec- 
tron temperature ti within the [0 11] zone in high-metallicity 
(12±log(0/H) > 8.25) Hn regions. The basic idea is the fol- 
lowing. If we apply the equation used to calculate the ionic 
abundance ++ /H + not to the entire Hll region but only to 
the O zone, then this would yield, not the ionic O++/H+ 
abundance but the total O/H oxygen abundance instead. 
We require that the equation for ++ /H + applied to the 
++ zone and the one for + /H + applied to the + zone 
result in exactly the same value of the oxygen abundance. 
This condition allows us to derive a relation between t2 and 
t3. We have applied this method to a sample of 372 Hll 
regions selected to have high-quality m easurements by us- 
ing the ff relation JPilvuein et al.ll2006l) . We find that the 
correlation between t2 and t3 is tight and can be approxi- 
mated by a linear expression. The so derived t2 - ts relation 
is independent of photoionization models of H 11 regions. 

The derived relation can be used to determine t2 and 
accurate abundances in high-metallicity H 11 regions with a 
measured t3 tempera ture. It can be also us ed in conjunction 
with the ff relation of lPilvuein et all (120061) for the determi- 
nation of the t3 and t2 temperatures and oxygen abundances 
in high-metallicity Hll regions where the [OIII]A4363 auro- 
ral line is not detected. 

Our t2 - t.3 re l ation is close to the widely used relation 
of lCampbell et aT] <ll986T) . However, the abundances derived 
with our t2 - t3 relation are slightly higher (up to ~ 0.1 dex) 
than those derived with the lCampbell et al] dl986l) relation. 
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